Non-linear feedback effects in coupled Boson-Fermion systems 
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We address ourselves to a class of systems composed of two coupled subsystems without any 
intra-subsystem interaction: itinerant Fermions and localized Bosons on a lattice. Switching on an 
interaction between the two subsystems leads to feedback effects which result in a rich dynamical 
structure in both of them. Such feedback features are studied on the basis of the flow equation 
technique - an infinite series of infinitesimal unitary transformations - which leads to a gradual 
elimination of the inter-subsystem interaction. As a result the two subsystems get decoupled but 
their renormalized kinetic energies become mutually dependent on each other. Choosing for the 
inter- subsystem interaction a charge exchange term, - the Boson — Fermion model - the initially 
localized Bosons acquire itinerancy through their dependence on the renormalized Fermion disper- 
sion. This latter evolves from a free particle dispersion into one showing a pseudogap structure 
near the chemical potential. Upon lowering the temperature both subsystems simultaneously enter 
a macroscopic coherent quantum state. The Bosons become superfluid, exhibiting a soundwave like 
dispersion while the Fermions develop a true gap in their dispersion. The essential physical features 
described by this technique are already contained in the renormalization of the kinetic terms in the 
respective Hamiltonians of the two subsystems. The extra interaction terms resulting in the process 
of iteration only strengthen this physics. We compare the results with previous calculations based 
on selfconsistent perturbative approaches. 

PACS numbers: 71.10.-w, 74.25.-q, 05.10.Cc 



I. INTRODUCTION 

A wide class of problems in solid state physics can be 
described in terms of interacting Boson - Fermion sys- 
tems. Examples are: 

i) the electron-phonon problem and its associated re- 
alization in superconductivity and polaron formation, 

ii) fermions interacting with spin fluctuations, relevant 
for the description of heavy fermion systems as well as 
for high temperature superconductors, 

iii) certain systems which can be mapped into Boson- 
Fermion interacting systems via Hubbard Stratanovich 
transformations, 

iv) the Anderson impurity and Kondo problem and fi- 
nally 

v) the Boson-Fermion model (BFM) for high T c su- 
perconductivity, believed to describe a coupled electron- 
phonon system in the cross-over regime between weak 
and strong coupling. 

In order to obtain the correct low energy physics in 
these various scenarios of Boson-Fermion interacting sys- 
tems the mutual feedback effects caused by the interac- 
tion between the two subsystems must be handled prop- 
erly. It invariably gives rise to effective time dependent 
interactions among the constituents in each subsystem 
which quite generally can be obtained by the standard 
field theoretical method based on functional integrals. 
This method is particularly suited if one considers the 



limit of infinite dimensions where it has been developed 
in great detail and is known as dynamical mean field 
theory p|. The mutual feedback effects have been re- 
cently studied on the basis of this method for the Boson- 
Fermion model 0] and for the many polaron problem 
in the Holstein model ||. While this method is non- 
perturbative and capable of describing the low energy 
physics, it is, as up to now, restricted to the study of 
local quantities. In situations where the dimensional- 
ity and anisotropy of the physical systems play a role, 
such as believed to be the case in the high T c cuprates, 
one has to resort to different techniques to handle these 
feed back effects. There are many physical situations 
in which already a relatively small inter-subsystem cou- 
pling leads to substantial fluctuations in each of them 
and which hence requires the selfconsistent determina- 
tion of the mutual non-linear feedback effects between 
the Fermions and Bosons. In those cases perturbative 
approaches in form of selfconsistent diagrammatic tech- 
niques can and have been applied for the electron-phonon 
problem (the Eliashberg approach B) and for the BFM 
H . These approaches however totally neglect vertex cor- 
rections. 

Controlled perturbative methods in the spirit of renor- 
malization group techniques have been recently pro- 
posed, the similarity renormalization scheme [pj and the 
flow equation technique [Q. These techniques are capa- 
ble in reformulating such interacting systems in terms 
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of renormalized Hamiltonians which capture the low en- 
ergy physics, which is achieved via an infinite series of 
infinitesimal unitary transformations. Contrary to the 
standard simple unitary transformations which treat the 
different energy scales in the problem in a single step (and 
therefore generally fail) these continuous transformations 
deal with each energy scale in a sequence of transforma- 
tions and by doing so are capable of extracting the low 
energy physics in which we are interested. 

This paper is organized in the following way. In the 
subsequent section II we review the essential points of 
the flow equation technique and apply it to the Boson- 
Fermion model. In section III we discuss the excita- 
tion spectrum of the renormalized Hamiltonian and com- 
pare the results to those previously obtained by different 
methods. In section IV we give a preliminary discussion 
of the superconducting phase properties of the BFM. 

Apart from testing this flow equation method for this 
model, this present study on the BFM goes beyond the 
studies so far reported in the past. For the first time we 
are able to treat on the same footing the normal and the 
superconducting phase and handle the question how the 
pseudogap evolves into a true gap below T c . 



II. THE FLOW EQUATION TECHNIQUE 

The method of continuous unitary transformation has 
been formulated by Glazek and Wilson || and indepen- 
dently by Wegner M . Instead of the single-step transfor- 
mation this method amounts to a procedure for simplify- 
ing the Hamiltonian's representation through the series 
of infinitesimal unitary transformations U(l), where I de- 
notes the continuous flow parameter. 

The continuous unitary transformation H(l) = 
U* (l)HU(l) gives rise to the following flow equation 



dHQ) 
dl 



= [r,(l),H(l)} 



(1) 



where r]{l) represents some arbitrary (anti-hermitean) 
generator 



(2) 



The choice of a specific form of rj(l) is usually de- 
termined by the physical situation under consideration 
In this paper we use Wegner's proposal JjJ 



r){l) = [H (l),H int (l)] 



(3) 



where H is the diagonal and H int the nondiagonal part 
of the Hamiltonian in a given representation. In gen- 
eral H int will be understood as the perturbation with re- 
spect to Hq. The generating operator Eq. (3) guarantees 
that under the continuous transformation the off diago- 
nal terms are monotonously reduced, eventually leading 



to a block diagonalization of the Hamiltonian, provided 
that no degenerancies are encountered Recently 
Mielke has proposed |8| some different form of the r)(l) 
operator which can be used for studying systems with 
degeneracies. 

One of the first condensed matter problems analyzed 
with use of the flow equation was the electron-phonon 
Hamiltonian |10|. It has been shown that eliminating 



the electron-phonon interaction induces the effective in- 
teractions between electrons which are attractive in the 
whole Brillouin zone. Near the Fermi surface this attrac- 
tion is strongly enhanced but it never becomes divergent 
as in the case for the classical Frohlich transformation. 
The method has also been successfully applied to a va- 
riety of other physical problems like the single impurity 
Anderson model |ll|] (the Schrieffer - Wolff transforma- 
tion has been improved), the strong coupling expansion 
for the Hubbard model Q (t/U expansion), the large 
spin Heisenberg Hamiltonian |13| (1/5 expansion), and 
for other topics such as dissipative quantum systems [fli) , 
light front QCD fl5||, quarkonium spectra p6[ and the 



Sine-Gordon model ||fj 

The main advantage of the continuous transformations 
is that an effective Hamiltonians can be derived which is 
valid not only in the low energy sector (like the standard 
Renormalization Group) but in the overall regime of en- 
ergies. In principle it is also possible to formulate the 
flow equations in such a way that lifetime effects can be 
studied too. Such type of flow equations have been used 
so far to account for dynamical effects in the context of 
spin-boson problem |f4j and also in the investigation of 
phonon damping effects due to the electron-phonon in- 
teraction Eal. 



A. Application to the Boson-Fermion model 

We shall now apply this flow equation technique to the 
BFM described by the following Hamiltonian 



H = H + Hi. 



(4) 



The free part (diagonal in the basis of the plane waves) 
consists of the kinetic terms of the Fermions and Bosons 

#o = YM - M^ka + J2( E « - 2 ^v ( 5 ) 

k,cr q 

The mutual interaction between both species is repre- 
sented by the charge exchange term 

Hint = ~7^J2 { v k,P b l+p c H c v] + v l. P b k+v c lAi) ■ ( 6 ) 



The flow equations control the evolution of the model 
parameters e G kl E q , Vk. q which get renormalized in the 
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course of the flow equations procedure. ^From now on 
we assume that they depend on the flow parameter I 
with the following initial conditions 



E q (l = 0) = A B , s%(l = 0) = e k , v klP (l = 0) 



(7) 



k,a q 

+ T7 E C/ P.fe.9(0 c pT c li c <?l c P+fe-?T + C (0 ( 9 ) 



p,k,q 



According to Wegner's definition (||) of the generating 
operator we have 



^(0 = ^E(4(0+4(0- wo 

x (^(O&fe+P^L - ^(O^I+pCfciCptJ (8) U p<k<g (l = 0) = , 



where the symbol : x :— x— < x > stands for the normal 
order product and c(l) denotes a c- number contribution 
to the Hamiltonian. We furthermore supplement the ini- 
tial conditions (]?]) with the constraints 



Upon iterating the flow equations (Jl]) new interaction 
terms are in general created which are not contained in 
the original Hamiltonian. Certain of those terms can be 
eliminated by reformulating the initial Hq in the follow- 
ing way 



AFD) 
l a,k 



(10) 

J2(E q -2»)4 BE K (11) 



k.a 



After some straight forward algebraic manipulations we 
finally obtain 



= "-iE **® ( v kA l ) b l+ P c kl^ + vtA l ) b k+p c tAl) 

k,p 

+ 4E ( a k,p( l ) + a q,k+p- q (l))v*k, p (l)v q .k+p~q(l)cl 1 cl l C ql C p+k . 



k,p,q 



8l 



k,p 



N 



k,p 



•M'JI n p+k '■ c k\ c k\ 



E [ a k,p( l )\ v k,p( l )\ 2 (- 1 + n { F k D) ) + a pM l )\ v pA l )\ 2n \ F k D) : b k+p b k+ 
k,p 

k,p 

+ M E b P+q b P+ k ( a kA l )+ a qA l )) V k,p( l ) V lA l ) C ll C ql + 
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N 



ABE) 
l fc+p 



k,p,q^k 



(12) 



where for brevity we introduce 

a k , p (l)=si(l) + sl(l)-E k+p (l). 
The expectation values are defined as 



< c\ a c ka >l 



i 



<&t 6g>i = n< BB >(l) = 



e (eJ(0-A«)/fcBT + 1 ' 
1 

e (E q (l)-2 fJ .)/k B T _ i 



(13) 

(14) 
(15) 



where rS FD \ n^ BE ^> are the Fermi-Dirac and Bose- 
Einstein distribution functions respectively. They depen- 
dent on I only through their arguments £&(/) and E q (l). 
In order to proceed with the numerical analysis of these 



a) our theory is valid up to the order v 2 , 

b) we omit any fluctuations of the form 

( c L C k* ~ ( C L C ^)) ( b l b P ~ ( b l b p)), 

c) as shown in the Appendix, contributions 
coming from the terms b P+q b p +kc\ a c qa for 
q =^ k arc of order 0(v 3 ). 

Within that procedure we finally arrive at the following 
set of the flow equations for the BFM 



flow equations (12) we shall neglect from now on the 
terms in the last two lines. Such a neglect implies that: 



dvk, P (l) 2 m m 

— j t — = -G<k,p( l ) v kA l ) 

de{{l) 



dl 



= ^E^p( / )i^p(oix+?(o 



(16) 
(17) 
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dejjl) 
dl 

dE k {l) 
dl 



dl 

rfc(Z) 
dl 



= ^E^(oi«m(oi 2 4+?(o 



2 

N 



E [afc_ Pi p(Z)|wfc- p ,p(0| 2 f- i + » 
p 

ap,fc-p|t' P ,fc-p(0| 2 «T F fc- ) p(0 



(18) 
(19) 



Let us concentrate from now on on exclusively two 
channels of the induced Fermion-Fermion interactions; 
the zero momentum BCS channel and the zero momen- 
tum density-density (d-d) channel respectively 



K, P (Z) + a 9 , fe+p _ 9 (/)] Vk iP (l)v q , k +p-q{l) , (20) 
2 



p.k 

v ■ 



(27) 
(28) 



^Eh^ (/)|u ^ (0|2 (" 1+ "S D)(0 ) 



(21) 



They denote specific elements U p B k CS ^ = U p ,- p ,- k (l — > 

oo ) and {/^ ''^ = U Pt k,kQ — * °°) which together with 
Eq. (|26|) become 



A formal solution for the flow equation ([16]) can be given 
right away 



(i) = v e -/ , [ e i( , ')+";c')-B*+p(i')] a «H' 



(22) 



but of course the dispersion functions and E p {l) 

ought to be determined selfconsistently via the other flow 
equations. 



B. Lowest order iterative solution 



The flow equation scheme is devised in such a way 
that the dominant renormalization takes place for the hy- 
bridization constant v k . p {l). To get some insight about its 
effectiveness we solve here the flow equations approxima- 
tively using on the right hand side of (|l6^ the bare 
(unrenormalized) energies £ k {l) — £k and E q (l) ~ A B 



{=E q =o). The resulting hybridization constant (E2|) re 
duces in this case to 



Vk,p(l) 



Vp,k 



(0 



-{e k +e p -A B ) 2 l 



(23) 



and has the desired property Vk, p (l — ► oo) = for all mo- 
menta, except when e k + £ p = A#. This situation cor- 
responds to a resonant scattering of two Fermions into a 
Boson state with their total energy being conserved. We 
shall see in section III that in the selfconsistent solution 
of the flow equations such a problem does not occur. Sub- 
stituting ( |23| ) into the flow equations (|l7| -pO|)and taking 
the limit I — > oo one obtains the renormalizcd quantities: 



(BE) 
4=0 



E { k R) = A B 



V 

I I N l 



AFD) 



(FD) 
k—p 



u 



(R) _ 



p,k,q 



-2A B 

(e fc + e p - A B ) 2 + (g g + £ fe+p _ 9 - A B ) 2 



(24) 



(25) 



(26) 



where el stands for the new Fermion spectrum, valid 
for both spins. 



U, 



(BCS) 



(Ai 



2e t 



2e k 



P j & 



(A B - 2e p ) 2 + (A B - 2e fc ) 2 



U. 



(d-d) 
p.k 



A, 



(29) 



(30) 



They are divergent in certain regions of the Brillouin 
zone, but as we shall show below, a selfconsistent nu- 
merical solution of the flow equations smoothens such 
divergences into a regular behavior. 



C. Comparison with standard perturbation the- 
ory 

Let us next discuss the results of the first iterative solu- 
tion of the flow equation method and compare it to that 
obtained by standard perturbative studies of the BFM, 
discussed in detail in Refs [§|,0. The second order ex- 
pansions for the Fermion and Boson selfenergies of the 
single particle propagators are given by 

T, F {k,iu) n ) = 

\v\ 2 k B T — ^ G F (p - k,iui m ~ iu> n )G B (p,iuj m ) 



1 



-\v\ 2 k B T j- G F {k - p, iuj m - iu n )G F (p, iuj n ) (31) 

P^n 

Our approximate solutions of the flow equations de- 
rived above (Section II B) were based on unrenormalized 
Fermion and Boson spectra and therefore are not self- 
consistent. Let us now determine the selfenergies using 
the bare Green's functions, as it has been done in Ref. 
H and which read 



£ F (fc,c) = M 2 ^E 
£ B (A;,w) = M 2 



AFD) 



ABE) 



n p-k + U q=0 

N^ui + h i-A B + e p - k ' 
1 1 _ n k FD p ] ~ n< p FD) 



N uj + 2/i — e p - k — s p 



(32) 



(33) 
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with n^o 3 = [exp{(A B - 2fx) /k B T} - l]" 1 . The effec- 
tive quasi-particle spectra are then given by the solutions 
of following equations 



III. NUMERICAL SOLUTION OF THE FLOW 
EQUATIONS 



2/i - A B 



E F (fc,c^)=Q 
S B (fc,wf ) = 



(34) 
(35) 



In the limit of small v we can put 



2/i 



and then substitute these quantities into the expressions 
for the selfenergies (||3,E33|) which results in 



2/x 



N ^ 



(FD) 
% p-k 



ABE) 
l q=0 



£k + £p-k — Aj 



1 1 N ^ i 



AFD) 



AFD) 
l k-p 



£fc- 



(36) 
(37) 



The difference between the Fermion spectra (|36|) de- 
rived from this perturbative approach and that derived 
from the flow equation approach can be remeded af- 
ter having realized that the standard perturbative study 
describes an effective Fermion spectrum while the flow 
equations method, reformulating the Boson-Fermion in- 
teraction, results in a) a renormalization of the Fermion 
energies £k —* ei^ and b) an appearance of the Fermion- 
Fermion interactions. Taking both effects into account 

when evaluating the final Fermion quasi-particle spec- 

( /?) 

trum, the lowest order corrections to e k is given by the 
Hartree term; i.e. 



AR) 



AR) 



k,p,p p 



(38) 



In this section we present the results of the selfconsis- 
tent numerical solution of the flow equations. In order to 
solve the differential equations, Eqs. (|l6| - |l|) we imple- 
ment the Runge Kutta method. The I dependent phys- 
ical quantities are determined iteratively. Starting from 
their initial values, Eqs. (|^, [lO| , [ri| ) we determine them in 
the following step according to x(l + SI) = x(l) + Six' (I), 
where x'(l) stands for the derivative of x with respect to 
the flow parameter I, as given by the corresponding flow 
equation. The increment which we use for this proce- 
dure is the following: SI — 0.01 for I < 5, SI — 0.1 for 
5 < I < 100, SI = 1.0 for 100 < Z < 1000, and finally 
51 = 10 for 1000 < I < 10000 (the flow parameter I is ex- 
pressed in units of the inverse square bandwidth D~ 2 ). 
The major renormalizations take place up to I ~ 100 or 
500, but certain parts of the Brillouin zone are slightly 
affected by a further increase of I up to few thousands. 
We control our choice for the upper limit of I by: a) look- 
ing on the asymptotic behaviour of all the renormalized 
quantities, and b) by checking whether all the hybridiza- 
tion matrix elements Vk,q(l) decreased below 0.1% of their 
initial value. 

In order to compare the results of this method with 
the results previously obtained by a selfconsistent per- 
turbative treatment || we choose the same set of model 
parameters as used in the above mentioned previous stud- 
ies, i.e., v = 0.1, A B = -0.6 and n tot = J2 a ( c L c i°-) + 
2 (b\bi^ = 1 and consider a one-dimensional tight bind- 
ing structure with the initial dispersion Ek = —2tcos(ka) 
(we set the lattice constant a = 1 and use the bandwidth 
D = At as a unit throughout this work). 



From the approximate solution (26) we have 



U, 



(R) 
k,p,p 



£k + £p- 
which eventually leads to 



AR) ^ AR) 



AFD) 



N ^ e k +e p ~ A E 
1 n (BE) + n {FD) 

i I |2 1 9=° p 

p * 



(39) 



(40) 




and which thus is identical to the expression, Eqn. ( pq ) 
obtained from the diagrammatic perturbation theory 
analysis. 



FIG. 1. Sketch of the hybridization matrix Vk, q (l) for 
I = 100. 
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It is instructive to first of all have a look how the hy- 
bridization matrix Vk, q (l) evolves in the course of this 
renormalization technique. In Fig. 1 we show Vk,q(l) 
for I = 100. Most of the terms of the matrix (k,q) are 
practically reduced to zero while there is a region in the 
momentum space for which the hybridization is only very 
weakly affected. 

To understand why a situation like that shown in Fig. 
1 takes place let us come back to the approximate solu- 
tion of the flow equations, without treating the Fermion 
and Boson spectra selfconsistently (section II B). By in- 
spection of Fig. ( p3| ) we see that for k, q such that 
£k + £q — As = the coupling constant Vk, q {l) remains 
unrenormalized. In order to determine around which mo- 
menta (q, k) this happens, we use Ek — ~2t + tk 2 . This 
gives us a topology of concentric circles around a mean 
radius given by 



k* = 



2a 1 



Ab 
D 



0.894427 . 



(41) 



Solving the flow equations selfconsistently we do not en- 
counter such a pathological situation but still, near those 
q and k points, the renormalization of Vk,q(l) evolves only 
very slowly. There are different characteristic lo points 
from which efficient renormalization starts for the mo- 
mentum sector near the resonant scattering Ek + e q = 
Ek+g (as also pointed out before by Ragawitz and Weg- 
ner 1 18 1 for the electron-phonon problem) and away from 
it. 

Fig. 2 shows the evolution of the hybridization con- 
stant along the q = — k cross-section as a function of 1. 
It is clear that renormalization of all the model param- 
eters is necessary, otherwise the total elimination of the 
Boson-Fermion interaction is difficult or even impossible 
to fullfill. 



To keep the total number of particles fixed we have to 
tune the chemical potential. There are two effects ob- 
served in the behavior of the chemical potential. First, 
with a decrease of temperature, the chemical potential 
approaches from below the bottom of the Boson band. 
Simultaneously, the bottom of the Boson band lowers 
and this is the reason why below a given temperature 
(T ~ 0.007) /i starts to be pulled down. Of course, the 
relative distance Eq — 2/i is a monotonously decreasing 
function of temperature. 

In this 1 — D case studied here, condensation of Bosons 
can of course not take place. The chemical potential ap- 
proaches asymptotically the lowest Boson level E q —o but 
never touches it except at T = (see Fig. 3). 








0.01 0.02 0.03 0.04 



FIG. 3. Variation of the chemical potential as function of 
temperature (at I = oo) for As = —0.6, n to t = 1. The Boson 
band is shown as a shaded area. 



B. The Boson spectrum 
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FIG. 2. Evolution of the momentum dependence of the hy- 
bridization constant V-k,k(l) as a function of the continuous 
iteration parameter I — 10 Q . 



A. The evolution of the chemical potential 



Due to the interaction with Fcrmions the initially local- 
ized Bosons acquire itinerancy. The effective interaction 
between Bosons, being of order t> 4 , is hence neglected 
here. In Figs. 4,5 we show how, with a decrease of tem- 
perature, the Boson band becomes broadened and the 
effective mass of the Bosons tub gets reduced, finally 



0) 



r /4. As a relative 



saturating around ms(T 
quantity we use m° F — which refers to the bare 

initial Fermion mass for the ID tight binding case. A 
parabolic curvature of the long wavelength limit k — ► 
is characterized by the inverse effective mass. 

For all temperatures studied by us the Boson disper- 
sion function exhibits a resonant like feature (a kink) 
around momentum q = 2k* ~ 1.8. This does not cor- 
respond to the value of 2kp as was initially mistakenly 
believed 0. By inspecting the location of such a kink 
for other sets of model parameters we conclude that it is 
mainly depending on the choice of As ■ This kink occurs 
for such momenta which satisfy the condition 



G 



(42) 



where q is a wavevector in the first Brillouin zone. Since 
within a precision of the order v 2 the value of the Boson 
energies is around the initial As, k* is roughly deter- 
mined by the condition 2ek* = A^. In an approximate 
study of the flow equations based on a first iterative sub- 
stitutions (section II. B) one can check that the function 
E q , Eq.(p5[) in fact diverges for q — 2k*. 



0.05 




n/2 k n 

FIG. 4. The effective Boson dispersion for a chosen set of 
temperatures. Notice the increase of the band width with 
a decrease of temperature. At very small temperatures the 
bottom of the Boson band asymptotically approaches the po- 
sition of 2/x and as a result the top of the Boson band gets 
somewhat pushed up. 



m B (T) 








0.01 



0.02 



T 



0.03 



FIG. 5. Low temperature behaviour of the effective Bosons 
mass in units of the initial Fermions mass m° F . In the insert we 
draw the long wavelength limit of the Boson spectrum for the 
three lowest temperatures evaluated; T = 0.001, 0.005, 0.007. 



ters we obtain an effective interaction between Fermions. 
Figs 6,7 below illustrate the momentum characteristics 
of the two channels defined in Eqs. (27,2§|). 

Again, we notice certain characteristic features appear- 
ing for the momenta corresponding to k* . The interac- 
tion Ujf d ' has a rather regular behaviour: for momenta 
such that k 2 +q 2 < k* 2 this interaction is attractive, while 



elsewhere it is repulsive. Around the region k 



= k* 



we observe a changeover, which looks quite singular. 



U, 



The interaction U^ B q CS ^ shows a similar behaviour as 

{d-d) 



but only along the cross-section q = k (see the 
top of Fig. H). Otherwise, the corresponding change- 
over between the attractive and the repulsive interaction 
regimes has a smooth character (see the bottom of Fig. 

I). 



I i (BCS) 
U k,q 




FIG. 6. The BCS interaction strength U { k B ° S) for 
T = 0.001. 



1 1 (d-d) 
u k,q 




C. Interactions between Fermions 



As a result of the renormalization of the model parame- for T = 0.001. 



FIG. 7. the density-density interaction strength Ujf" q ^ 
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Temperature has a negligible effect on the effective 



Fermion-Fermion interaction. For example, U ( 



(BCS) 
0,0 



de- 



creases only by about 0.4% when the temperature is var- 
ied from 0.1 to 0.001. No qualitative change is observed 
at all. Nevertheless temperature is an important factor in 
as far as the effectiveness of this Fermion-Fermion inter- 
action is concerned, since with a decrease of temperature 
kp — > k* . To see that we refer the reader to consult Fig. 
H keeping in mind that £o - Ag. As seen from the Figs. 
8 the interactions for \k\ < k* have attractive character 
and are strongly enhanced (at least the elements Uk, q =k) 
infinitesimally below the momentum k* . One would nat- 
urally expect strong effects of these interactions if the 
Fermi vector was located just below k* . 



u 
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FIG. 

(d-d) 



8. Cross-sections of the interactions Ul B q CS ^ and 



k along a) q = k when both interactions are identical 
(top figure), and b) for q = (bottom figure). In general 
both potentials are repulsive for momenta k > k* and be- 
come attractive if k < k* . 



D. The Fermion spectrum 



Studying the Fermion spectrum is a rather complicated 
issue because on one hand it is affected directly through 
the flow equations (17,18) and, on the other hand, indi- 
rectly through the induced interactions as discussed in 
the preceding section. The effect of temperature turns 
out to be a very important factor in this problem. In the 
physical regime of interest (i.e. when the effect of the 
Boson - Fermion coupling manifests itself strongly) the 
Fermi momentum is situated just below k* . Upon de- 
creasing the temperature kp moves closer and closer to 
that value. As a result Fermion - Fermion interactions be- 
come more and more effective and eventually may lead to 
destruction of the quasi-particle nature of those Fermions 
in this model. 

Below, we discuss what kind of fermion spectrum arises 
purely on the basis of the renormalization of the free par- 
ticle energies given in Eqs. (17,18). In Fig. 9 we plot 
the dispersion ej^ versus k and in the inset its deriva- 
tive with respect to initial Eh- Using this derivative we 
compute the density of states (DOS) 



p(«) = ^E*(«+/*-^ 

A,- 

= jdep\e) 



k J 
l 



(43) 



where p°(e) — N 
of the Fermions. 



1 *C/c 5{e — Ek) is the initial bare DOS 




-0.5 







TU/2 



FIG. 9. The renormalized (I — > oo) Fermion spectrum e^. 
for low temperature T = 0.005. Notice the tendency to form 
a gap around the Fermi vector. 

We notice that below a characteristic temperature T F 
(at which the chemical potential starts to be pinned at 
the bottom of the Boson band) there appears a pseu- 
dogap centered around the Fermi level. In our case 
T F ~ 0.125 in units of the initial bare Fermion band- 
width. 

An important question is the qualitative change-over 
of the pseudogap into a true BCS type gap below the con- 
densation temperature T c . This question can be tackled 
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in a reasonably controlled way in the ID case. T c is then 
identically zero and at T = all Bosons are condensed 
in the q = state. Remember that we are dealing here 
with an effectively free Bose gas on a lattice, the Boson- 



Boson interaction 
thus obtain 



of order v - being neglected. We 



de k {l) 
dl 

-k{l) 



4n?(0 ( £fc (Z)-M(0)K-fc(0| S 



dpll) 
dl 



dl 

1 dE (l) 
2 



= -4(e fc (l)-/i(0) «*,-*(0 , 
£fe(0 - MOI K,-fc 



1 



which follow from the general flow equations ( |l6[|l9| ) in 
the T = limit. tiq{V) denotes the concentration of con- 
densed Bosons and /i(Z) = E (l)/2 in this limiting case. 
We need not find the whole Boson spectrum E q (l) be- 
cause finite momentum states E q ^Q are not relevant in 
the ground state (at least in absence of some external 
fields). 




-0.2 0.2 0.4 0.6 0.8 CO 

FIG. 10. The DOS for the renormalized Fermion spectrum 
@ for T = 0.02 (dotted line), T = 0.01 (dashed line) and 
T = 0.005 (solid line). Upon lowering the temperature a pseu- 
dogap forms near the Fermi energy i.e., (uj = 0), caused by 
the bonding two- Fermion states [22]. Some modification of 
the DOS is also seen for higher energies which corresponds 
to the antibonding two- Fermion state [22]. In the inset we 
compare the results obtained within the present flow equa- 
tion technique (solid line) with those derived by selfconsistent 
diagrammatic techniques [5] (dotted line) for T = 0.005. 

The set of Eqs.Q,|^,|6]) is rather straightforward 
to study because the momentum dependence of the 
hybridization constant Vk,-k(l) enters only through 
Hence one can use the effective DOS p(e,l) = 
l/NJ2k^( £ ~ £ fe(0) m Eq. (Eq). We solve numerically 
these equations for a fixed total concentration of particles 
which is subject to the supplementary condition 



«o S (0 



-ntot 



de(l) p{e,l). 



(47) 



(44) 
(45) 



I 
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2 (0 (46) 



=L -0.02 
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S 0.21 

O 

C 0.20 



500 / 1000 

FIG. 11. Evolution of the renormalized chemical potential 
(top) and concentration of Bosons Uq (I) (bottom). Total 
concentration is kept fixed at ntot = 1. 
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-0.5 
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FIG. 12. Difference between the DOS of the interacting 
p(uj) and of the free system p°(co) for n to t = 0.6597 and tem- 
peratures displayed in the legend. 

Fig. [ll] shows how the chemical potential is renormal- 
ized through this continuous transformation technique. 
Notice that its variation is of the order v 2 , as in the case 
of finite temperatures. The total concentration of Bosons 
(which are all in the condensate state at T = 0) is rather 
weakly affected by renormalization (of the order of 4 %). 

The important outcome of this T — case analysis 
is seen in the spectrum of Fermions. The asymptotical 
solution at I — > oo yields a true gap in e)^ which is 
formed around the chemical potential. The size of this 
gap is in very good agreement with the mean field theory 



prediction, i.e. A(T = 0) 



(see Fig. h3h 



In Fig. |12j we summarize our results obtained sofar in 
this section and which permit us to make some conjec- 
tures as to the evolution of the pseudogap into the true 
superconducting gap as the temperature is lowered. First 
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of all we notice two distinct energy scales which define 
these two gaps: 

i) the superconducting gap being very sharp and being 
controlled by the first power in the coupling constant v 

ii) the pseudogap, being evident in form of humps 
whose positions slightly move closer together as the 
temperature decreases, varies with the second power of 
the coupling constant v (remember that in the normal 
state our flow equation approach the results obtained are 
within the order of v 2 ). 

The relative size of the pseudogap and the zero tem- 
perature true gap can vary considerably as a function of 
the Boson concentration, as can be seen from Fig. |l3|. In 
Fig. |lj we have chosen a situation close to the experi- 
mental situation in the high T c cuprates where the two 
gaps are of comparable size. 

In the next section we shall show that this same tech- 
nique, applied to the superconducting state, gives the 
superconducting gap varying with v. It is clear that this 
difference in the size of these two gaps does not mean 
that they are of different origin! We shall come back to 
this question in the last section of this paper. 



0.15 



Q. 




0.05 



FIG. 13. Variation of the true gap at T = as function of 
Boson concentration obtained by the flow equation method 
(solid line) compared with the mean field solution (dashed 
line). The filled circles show the corresponding magnitude of 
the pseudogap obtained by the flow equation method for a 
normal phase at T = 0.01. 



IV. THE SUPERCONDUCTING PHASE 

Bosons are not able to condense at any finite (nonzero) 
temperature unless the dimensionality of the system is 
higher than two. In this section we try to reach some 
preliminary conclusions regarding the superconducting 
phase of the BFM on the basis of flow equation method. 
We assume that at least some fraction of the Bose subsys- 
tem is in the condensed state; in other words we consider 
T < Tbe- We show that the appearance of such a con- 



densate is inevitably related to the formation of a true 
gap in the Fermion subsystem. Our estimation of this 
gap is in agreement with the mean field theory result for 
this model pjj . 

A. The Fermion Subsystem 

Given the existence of a certain fraction of condensed 
Bosons we put the chemical potential at the q — level 
of the Boson energy spectrum, i.e., = Eg =0 {l)/2. In 
the flow equation for Fermion energies, Eqs. (O) we then 
have a dominating contribution coming from the n q=0 
thermal factor and, as a consequence, can simplify this 
equation to (which for T = is exact) 



de k (l) 
dl 



4n B (T,l)(e k (l)-n(l))\v k ^ k (l)\ 2 



(48) 



n B (T,l) denotes the concentration of condensed Bosons 
at temperature T. This equation is coupled to the 
flow equation for the hybridization coupling v k ^_ k , given 
in Eq. (|45|) . By inspecting the bottom Fig. [n] we 
see that the I dependence of the condensate concen- 
tration can be dropped, n B (T,l) ~ n B (T). For the 
strictly three-dimensional system the temperature evo- 
lution is given through the standard relation n B = 

n B 1 — (T/Tbe) 3 ^ 2 , where n B is the total concentra- 
tion of Bosons. Similarly, the / dependence of chemical 
potential can be dropped because of the following argu- 
ments: 1) for momenta close to the Fermi surface, the 
renormalizations of e k (l) are of the order of vy n B (T) 
while the chemical potential undergoes a renormaliza- 
tions of the order v 2 (see top Fig. [ll])) for momenta 
which are far from kp (i.e. in the high energy sector, 
using a terminology of the renormalization group the- 
ory) renormalization is not effective at all (except for the 
nevralgic k* point which is irrelvant for the purpose of 
the present discussion). 

Thus, without loss of generality and precision, we can 
rewrite the flow equations in the following from 



dl 
dv(l) 
dl 



-4£ 2 (ZM0 , 



from which immediately follows 

n B v 2 (l) +£ 2 (0 = const. 



(49) 
(50) 

(51) 



where t;(l) — e(l) — u measures the Fermion energy from 
the chemical potential. It is evident from Eq. (50) that 
for any nonzero £ the hybridization must evolve to zero 
in the infinite / limit, i>(£ ^ 0) — > 0. Consequently the 
renormalized spectrum becomes 



,(R) 



fi = sign (e 



m) 



2 . B 2 

+ rift v z 



(52) 
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The gap formed around the chemical potential in the 
Fermion subsystem is given by 



,(T) = v^rj 



(53) 



which is the same as predicted by our previous stud- 
ies Jl9|,^l| of this model. This result confirms that the 
BFM is characterized by a single transition temperature 
T c at which the Bosons start to condense and Fermions 
to form a gap in their excitation spectrum. Hence 



rp rp{B) rp 

1 c — 1 BE ~ 1 > 



BCS- 



B. The Boson Subsystem 

A very important issue in this context is to under- 
stand the impact of the gap in the Fermion spectrum on 
the excitation spectrum of the Bosons. Unfortunately 
the general flow equation, Eq. (|l^) can not be handled 
analytically, not even in the small q limit. 

In order to get some insight we study numerically this 
equation together with the constraint, Eq. ( |5(i| ) which, 
evidently, applies only for T < Tbe and for dimensions 
larger than two. The set of equations we have then to 
solve are Eqs. (19,22) together with 

MO - M (Z)] 2 = [e k - m] 2 + n$ [v 2 - k,_ fe (Z)| 2 ] (54) 

which is a direct consequence of the constraint, Eq 
(51). We fix the chemical potential at the level fi(l) = 
E q=0 (l)/2 throughout the iterative solution procedure. 
Eq. (Q) controls the formation of the gap in the Fermion 
spectrum as a result of the presence of Boson in the con- 
densate. With the constraint, Eq. (|54|) included, we in- 
vestigate the Boson spectrum E q by restricting ourselves 
to performing a one-dimensional momentum summation 
in Eq. (19). Such an approximative procedure reduces 
considerably the numerical complexity and is expected to 
give at least qualitatively correct results for dimensional- 
ity larger than two. A complete numerical study on this 
point will be reported in some future work. 




FIG. 14. Momentum characteristics of the Boson exci- 
tation spectrum in the superfluid state of the BFM with 
ntot = 1. Solid lines correspond to various temperatures as 
marked and the dotted line is taken from the selfconsistent 
solution for the normal phase at T = 0.005. Inset: the same 
function for the superfluid (solid line) and normal phase (dot- 
ted line) within the half Brillouin zone at T = 0.005. 
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FIG. 15. Sound velocity and the condensate concentration 
versus T for the superfluid phase of the BFM with n to t = 1. 



/From the numerical analysis of the above flow equa- 
tions we obtain E q . In the long wavelength limit 
(q — > 0) it significantly deviates from its behavior, ob- 
tained above, in the normal state. In order to illustrate 
this, we plot in Fig. [li] the momentum dependence of 



E, 



(R) 



Eq=Q 



0.1 0.2 0.3 0.4 0.5 q 



j q for several temperatures. Clearly the 
curves show linear behavior up to momenta q ~ 0.15 
and, what is more important, show a non-zero crossing 
point with the ordinate. Its value determines the sound 
velocity v s (T) and marks the presence of a collective ex- 
citation in the superfluid Bose subsystem. That such 
a sound wave mode is completely absent in the normal 
phase can be seen by the dashed line in Fig. [l4| which 
crosses the ordinate at zero. 

By inspection of Fig. |TJ one notices the decrease of the 
sound velocity with increasing temperature. Simultane- 
ously the region of the q 2 behavior of the spectrum starts 
to shrink. In Figs. (15,16) we show the dependence of the 
sound velocity versus temperature (for a total concentra- 
tion of carriers n tot — 1) and versus total concentration 
at T = 0.0. These results agree well with the predictions 
for the BFM obtained earlier by means of the dielectric 
function formalism In the so-called Bose limit, i.e. 
when the concentration of Bosons is not small, the sound 
velocity has been shown to gradually decrease with an 
increase of temperature towards T c . On the other hand, 
the sound velocity of the ground state gets reduced when 
the total concentration of carriers increases (see Fig. Il5] 
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of our present calculations and compare with Fig. 2 of 
Ref. pl|). Only in the dilute regime for small Boson 
concentrations (the so called BCS limit) one expects a 
behavior, qualitatively different from that studied here. 



0.1 







^n =0.01 

i 




i 
i 

i°- 05 

^ 0.25 0.5 


1.0 



FIG. 16. Sound velocity of the ground state T = as a 
function of total concentration n to t for Ab ~ —0.6. Notice 
that critical concentration above which superfluidity (super- 
conductivity) can arise is roughly equal n cr ~ 0.49. 

V. CONCLUSION 

Our previous studies j|] of the BFM indicated that the 
superconducting features of it arise due to the initially 
localized Bosons becoming itinerant, an effect which is 
triggered by a mutual feedback effect between the two 
subsystems. Concomitantely with this occurs the forma- 
tion of a pseudogap in the DOS of the Fermions, as the 
temperature is lowered and which in turn permits the 
Bosons to acquire longer and longer lifetimes. The rea- 
son for that is a reduction of scattering processes due to a 
diminishing number of Fermionic states available in this 
energy region. The opening of the pseudogap in those de- 
scriptions is linked to a renormalized Fermion dispersion 
which becomes fiat as the Fermi energy is approached 
from below, but, with at the same time, a substantial 
loss in spectral weight and lifetime broadening J24J . The 
end result is in effect a separation into two separate sub- 
systems, the Fermionic and the Bosonic one with their 
proper dynamics. Yet, due to the exchange coupling, 
they are mutually dependent on each other which leads 
to a single critical temperature, describing the onset of 
superconductivity, in both of the two effective uncoupled 
subsystems. 

The flow-equation technique studied in this work, 
renormalizes this inter-subsystem coupling to zero and 
hence is capable to make this interdependence of the dy- 
namics of these two subsystems explicit. The various 
results obtained here are correct to second order in the 



initial unrenormalized inter-subsystem coupling constant 
v. The effective Bosonic subsystem behaves essentially as 
a lattice gas of free Bosons with a temperature dependent 
mass. The effective Fermion subsystem shows a disper- 
sion which, upon approaching the Fermi wavevector rises 
almost vertically to some value above the Fermi energy 
within a very small regime around kp. This reflects a 
pseudogap structure in the DOS which is of order v 2 . In 
this regime of energies the effective renormalized intra- 
subsystem interaction is very singular (see Figs. (6 -@)) 
and is expected to give rise to the lifetime effects and the 
reduction in spectral weight seen on our previous stud- 
ies ||. For wavevectors greater than kp the effective 
Fermion dispersion remains quasi unrenormalized. 

The use of the flow equation technique enables us to 
treat the superconducting and the normal state on the 
same footing and, on the basis of that, to make some 
conjectures as to how the pseudogap evolves into the 
true gap in the superconducting phase. For a ID sys- 
tem where the superconducting phase is realized at T = 
the pseudogap evolves in from of a U-shaped curve which 
deepens until it touches the zero density level upon de- 
creasing the temperature towards T = + . Upon en- 
tering the superconducting state at exactly T = 0, this 
V shaped curve changes abruptly into a more conven- 
tional U shaped curve, known from standard BCS type 
superconductors. The pseudogap is characterized by two 
distinct humplike features in the DOS whose positions 
get slowly closer to each other as the temperature is de- 
creased. The energy difference between those two humps 
being of the order of v 2 . These pseudogap features are 
distinctively different from the gap structure in the su- 
perconducting state, as can be seen from Fig. 12. The 
superconducting gap shows a different variation with the 
inter-subsytem coupling constant, varying as v. The sizes 
for the pseudo- and the superconducting gap can in prin- 
ciple be quite different, as shown in Fig. |l3|. This does 
not mean that they are of different physical origin. We 
know from our previous studies 0] that the pseudogap 
is very much independent on dimensionality, while the 
superconducting gap evidently is dependent on it. The 
present study, further, suggests that the variation of the 
two gaps with the concentration of Bosons varies in op- 
posite direction. This can lead to a situation of a coex- 
istence of two gap like structures in the superconduct- 
ing phase, i.e., a superconducting gap and a remnant 
of the pseudogap. To what extent we should consider 
the pseudogap as a precursor of the superconducting is a 
question of semantics. For the model system considered 
here, as well as for the real High T c cuprate materials 
the pseudogap is caused by amplitude rather than phase 
fluctuations. Amplitude fluctuations are a prerequisite 
of the superconducting state but are not by themselves 
sufficient to guarantee its materialization. As a conse- 
quence, upon reducing the dimensionality of the system, 
the superconducting state is suppressed and tends to an 



12 



insulating state with a characteristic upturn of the resis- 
tivity at low temperatures @. 

The study of the change-over between the pseudo- 
gap and the superconducting gap for the more realistic 
anisotropic 3D case, together with a careful study of the 
lifetime effects controlled by the intra-Fermion subsys- 
tem interactions is presently under investigation and will 
be reported in some future study. 

Finally, for reasons of completeness, we should men- 
tion the theoretical studies of the pseudogap phenomenon 
based on the so-called BCS-Bosc Einstein cross-over sce- 
nario. Such studies are based on effective BCS type 
coupling as well as the negative U Hubbard model p5| . 
Within this approach similar questions of a change-over 
form the pseudogap into a true gap have been considered 
p6|| . A discussion the precise differences between this 
scenario and the BFM scenario as concerns the physics 
and applicability to the High T c cuprates lies outside the 
frame of subject discussed in the present paper. 
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Appendix 

Let us consider the following modification of the gen- 
erating operator 



(2) 



(55) 



where rj is given by equation (pj) and where we choose 



.,(2) 



E (lp.k, q (l)bl +q b P+k cl a c qa -h. 



<y p,k 7 q^k 



(56) 



The coefficients j° k Al) can be selected in any arbitrary 
way provided that Vk, P (l — > 00) — > still holds. 
By a straightforward calculations one verifies that 



(7 p,k,q^£k 

J 



P +q U P+k L k<j C 1 cr 



cLc 0CT + h.c 



O^v^U). (57) 
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If we now use the Ansatz 

e\ 

fi) = 



+ el - E, 



p+k 



E 



p+k 



E 



>p,k,q 



(0 



.T 



E, 



p+k 



E,, 



E r , 



p+k ^p+q 

then we effectively obtain 
1 

N 



i u k,p u q,p 



I v p,k v p,q 



(2) 



H 



E 

p.k,q^k 



bl+qbp+k 



[(«fc,P 

X 4j. C ?i + ( a P,k + a P,q) V p,k V P,q C kl C <lt 



q,p) v k,p v q,p 



(58) 



(59) 



'V, 



The terms containing become of the order 0(v 3 ) and 
from the flow equation ( |20| ) we can estimate U ~ w 2 
which yields that 7?/ ~ 0(i> 4 ). 

Thus on the right hand side of equation (^) we ob- 
tain the same term as in ( |l2| ) but with an opposite sign. 
These terms subtract each other if one uses fj in the flow 
equation instead of the initial one (^) . The modified con- 
tinuous unitary transformation does not generate any in- 
teraction of the form bp +q b p +kc\. a c qa for q 7^ k unless 
hybridization constant v is large enough. 



0{v 3 ) . (60) 
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